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Previous results on the statistics of avalanche detectors are generalized to 
the case where electrons and holes suffer collision ionizations with unequal 
probability. It is assumed here that the ratio of collision ionization prob- 
abilities per unit length of weaker-to-stronger carrier is a constant k in- 
dependent of position in the high-field region. The moment-generating 
function of the random avalanche gain G is obtained as a function of k 
and the average gain G, and is used to obtain Chernov bounds on error rates 
of digital optical receivers employing avalanche detectors. It is shown the 
required energy per pulse to achieve a given error rate decreases as k de- 
creases for fixed G. For each k > 0, there is an optimal mean gain G op , 
resulting in minimum required energy per pulse. At k = 0.1, G opt ~ 100 
and the required energy is within 10 dB of that required with very high 
gains (a few thousa?ul) at k = 0. 

I. INTRODUCTION 

In a previous paper 1 results on the statistics of two particular ava- 
lanche detectors with applications to optical communication were 
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presented. It was required either that only one carrier suffer collision 
ionizations in the high-field region (unilateral gain) or that both carriers 
suffer collision ionizations with equal probability per unit length in the 
high-field region. The present work allows for more general unequal 
ionization probabilities per unit length with the requirement only that 
the ratio of the two quantities be constant throughout the high-field 
region. The moment-generating function of the random gain is obtained 
as a function of this ratio and the average gain. The results are con- 
sistent with unpublished conjectures of R. J. Mclntyre. 2 The moment- 
generating functions are used to obtain Chernov bounds on the error 
rates of digital optical receivers employing avalanche detectors and 
using either coherent or incoherent light. Results on avalanche statistics 
are summarized in Section V. Numerical results on the Chernov bounds 
are given in Section VI (6.5) and Section VII. 

II. MODEL OF THE AVALANCHE DETECTOR 

The avalanche detector is a device in which thermally or optically 
generated hole-electron pairs generate additional hole-electron pairs 
through collision ionizations. Within the device there is a "high-field 
region" where holes have probability /3(x) per unit length (which depends 
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Fig. 1 — Avalanche detector. 
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upon the position x) of suffering a collision ionization as they travel to 
the left under the influence of the electron field (see Fig. 1). Electrons 
traveling to the right have a probability a(x) per unit length of collision 
ionization. Carriers can be created within the high-field region due to 
thermal effects or due to the presence of incident light. Carriers can 
also drift into the region if they are generated outside of the region. It is 
assumed that all collision ionizations are independent. This requires 
that the mean distance between ionizing collisions be large compared 
to the distance over which a carrier can randomize its momentum after 
a collision. 3 Hole-electron pairs created through collision ionization can 
in turn generate additional pairs by the same mechanism. This is the 
avalanche process. 

III. THE STATISTICS 

We seek the statistics of the random total number of hole-electron 
pairs which result ultimately through collision ionizations when an 
initial hole-electron pair is injected into the high field at some position x. 
Define p u (n, x) as the probability that n pairs ultimately result in- 
cluding the initially injected pair. The moment-generating function of 
the number of pairs M„(s) is therefore 4 

M (s, x) = ^p g {n,x)e n . (1) 

i 

We shall derive M (s, j,-). Before proceeding we must review some well 
known results which will be needed in that derivation. 

If ( Xt } are random variables which are independent, then the moment- 
generating function of the sum of the { a:,- j is the product of the individual 
moment-generating functions. 4 The semi-invariant moment-generating 
function SIMGF of a random variable X having probability density 
p x (x) is the natural logarithm of the moment-generating function of X 

Ms) - In [M x (s)] = In [ f) p x (x)e"j • (2) 

The SIMGF of a sum of n independent random variables is the sum of 
the individual SIMGF's. 

We can now proceed to derive M„(s, x). Divide the high-field region 
into K intervals of width dX = W/K. See Fig. 2. Label these intervals 
1, 2, 3, • • • j, • • • iv. If a hole-electron pair is injected into interval j, 
define, as above, the probability density of the total number of pairs 
ultimately resulting in the avalanche process (including the initial pair) 



3078 



THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1971 



as p g (n, x) where x is taken as the center of interval j. The hole of the 
initial pair moves to the left and the electron to the right toward 
x = or x = W respectively. As they pass through their respective 
intervals, new pairs may be created in each through collision ionization. 
We shall assume that the interval width dX is sufficiently narrow so 
that the initial pair carriers create either one or no new pairs in each 
interval. If the initial hole or electron generates a new pair in some 
interval k, then that new pair will ultimately generate N k pairs including 
itself through the avalanche process. Thus with each interval we can 
associate a number of pairs N k . This number equals zero if the appro- 
priate initial pair carrier suffers no collision ionizations in interval k. 
This number equals one or more if the appropriate initial pair carrier 
suffers a collision ionization in interval k. The total number of pairs 
ultimately generated through the avalanche process including the 
initial pair is one plus the sum of the {N k }. Since collision ionizations 
are all independent, all the N k are independent. Thus we have the 
SIMGF of the total number of pairs given by 



W/dX 

t„(s, x) = s + 2 fai') 



(3) 



where s is the SIMGF of the deterministic initial pair and \f/ Nk (s) is 
the SIMGF of N k . 

The SIMGF of N k is obtained as follows. The probability that 
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Fig. 2 — Avalanche process. 
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N k = is 1 — 7 dX where y = a(x) if interval k is to the right of interval 
j where the initial pair enters, 7 = (3(x) if k is to the left of ;. (x is eval- 
uated at the center of interval k.) The probability that N k = Z > is 
the probability of a collision ionization in interval k times the probability 
that a new pair created at interval k ultimately results in Z pairs 
including itself. That is 



Pr (N k = Z) = 7 dX p g (Z, x) for Z > 
with x evaluated at the center of interval k. Thus 



(4) 



ftr»(«) = In 



(5) 



where 



£ Pr (Y, = zy 

z = a 

= In [(1 - 7 dX)e'° + £ 7 dX p,(Z, x)e' z ] 

= In [I - ydX + y dX M„(s, x)] 
= ln[l - y dX + y dX e*' l,,x) ] 

7 = a(x) if k > j 

7 = 0(a:) if k < j 

x = center of interval k. 

Using (3) and (5) and taking the limit as dX gets infinitely small, one 
obtains 

+ I a(a-')(e^ (SI,) - 1) (to'. (6) 

Equation (6) is the critical equation for determining \f/ a (s, x), and thus 
il/ 5 (s, x) = exp [^„(s, a;)]. Using Leibnitz's rule for differentiation of 
integrals one obtains 



£*.(•,*) - DSfr) - a(x)](^' { - rt - 1). 



The solution of (7) is 
yp B (s, x) = In 



1 



1 - C exp I / \p(x') - a(x')] dx') 



(7) 



(8) 
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where C = (e*» <,,0) - l)/e*' <,,0) which can be checked by substitution. 
Substituting (8) into (6) one obtains the particular result 



Us, 0) 



= s + 



•'o 



1 -1 



dx'. (9) 



1 - C exp ( j" &3(a;") ~ «(*")] <&") 

If one makes the assumption 3 that at each point in the high-field 
region 

p(x) = k-a(x) (10) 

where k is a constant, one can solve (9) to obtain 

fcfa 0) = 8 - 5 + r^ In [ 



* - 1 LM (s, 0) - e tk - l) '[M.(8, 0) - 1]J 

= s + r ^ In [M„(s, 0) - e (i - l,J [M ; ( S) 0) - 1]] (11) 

where 

M B (s, 0) = e*' (,,0) and 5 = / a(x) da;. 

Jo 

One can write (11) in another way by making a substitution.* Define 
5'(s) implicitly by 

e s ' l,) = M B (s,0)e s -', (12) 

that is, 



*.<•.»> _ m /-o m = *—«+«'<•> 



= M e (s,0) = e'- ,+, ' (,) . (13) 



Using (12) in (11) one obtains the implicit equation 



--«'(.) = g.r e -« _ e -lj (14) 



which determines 5'(s) and thus M B (s, 0) through (13). 

Equation (14) is still not explicit. A numerical technique for solution 
is discussed in the next section. One can use (11) [or (14)] and (8) to 
obtain ^„(s, x) or M B (s, x) for any x. Recall that x is the point of entry 
of the initial pair. In the applications we shall be concerned with x = 
or x = W. That is, pairs are generated in a drift region outside the high- 
field region with carriers drifting into the high-field region. 

* Equation (14) will follow from (11) and (12) by tedious algebra. Further, (14) 
will not be used in the following results except to compare with Mclntyre's work 
in Appendix A. 
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The equation (14) is consistent with some unpublished conjectures of 
Mclntyre given in Appendix A. 

IV. NUMERICAL SOLUTIONS 

Equations (11) through(14) can be solved numerically. One technique 
is to differentiate (11) to obtain the result 

f s M.(8, 0) = M (s, 0)(^ J: )[l - I (M,(s, tye-V)*- 1 ] (15) 

where M B (0, 0) = 1. 

Equation (15) can be integrated with a computer to obtain M„{s, 0) 
explicitly. 

V. SUMMARY OF ANALYTIC RESULTS ON AVALANCHE STATISTICS 

From Sections I through IV and a previous paper by this author, 1 
we obtained the following: 

Assumptions: 

Holes travel toward x = 0, electrons toward x = W. 
Hole ionization probability per unit length = 0(x). 
Electron ionization probability per unit length = a(x). 
|8(x) = k-a(x), k a constant for all x. 
High-field region width = W. 

Definitions: 

c w 

8 = I a(x) dx, 

Jo 

p„(n, x) = probability that if an initial pair enters the high-field region 
at point x, n pairs will ultimately result through the avalanche process 
including the initial pair, 

G(x) = mean avalanche gain = ^ n Po( n > x )> 

i 

M B (s, x) = moment-generating function of p„(n, x) 

CO 

= H P (n, x)e' n . 

n-l 

Results: 
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i - *•<*• * = b- {M uS,o) 1} «p [<* - » />'> *F 

for all fc. (16) 
For /c = (Unilateral Gain) 

2a. M,(«, 0) = [1 - e 5 [l - e--'}]- 1 where (7(0) = e 6 . (17) 

2b. p,(», 0=| (^r 1 )" where G = (7(0). (18) 

For k ?* 0, k ^ 1 

3. £m„(s, 0) - Af,(«, 0)(^ lI )[l - I (M,(s, (Be""/)*- 1 ] ', 

where 

M„(0, 0) = 1 (19) 

where 

^) = (^T J )[ 1 -|^""]"' 

For fc = 1 (Equal Ionization) 

M„(0, 0) = 1 
where 

G = (7(0). (20) 

VI. APPLICATIONS TO RECEIVERS USING AVALANCHE DETECTORS 

6.1 General Comments 

We shall next apply the results of Section V to obtain bounds on the 
error rates of digital receivers. The receivers to be discussed here are 
the single- and twin-channel systems described below. We shall upper- 
bound the power required at the receiver to obtain a desired error rate 
using the Chernov bounds. 

6.2 The Receivers 

The twin-channel receiver is shown in Fig. 3. Depending on the state 
of a binary information source, one of two channels has optical output 
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power p for duration T while the other has output power p • EXT for 
duration T. EXT signifies extinction. EXT would ideally be zero, but 
is left finite and less than unity for practical reasons. The optical power 
falling on an avalanche photo diode causes the emission of photo- 
electrons which are multiplied along with the detector "dark current" 
through the avalanche gain mechanism. The detector outputs (multi- 
plied counts) are integrated in devices having thermal noises referred to 
their respective inputs. The integrator outputs are subtracted and the 
difference X is compared to a threshold of zero to decide what the 
information state was. 

The single-channel system is essentially half the twin-channel system 
as shown in Fig. 4. The single integrator output X is compared to a 
threshold y to decide upon the information state. 

6.3 The Chernov Bounds 

The Chernov bound is a useful tool for bounding the probability that 
a random variable X will lie above or below a given threshold 7. It is 
given as follows 5 

Pr(X> 7 ) &e*' M -'*' M | w ; ( ., 

provided s > 0, 
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provided s < 0, (21) 



where 



and 



fcr(s) = SIMGF of X = In IT p x (x)e' z dx\ 



*&) = f 8 M*)- 



6.4 Chernov Bounds for the Tioo Syste7ns 

6.4.1 Preliminaries 

We wish to determine the required power p to achieve a desired 
decision error probability for each of the receivers discussed in Section 
6.2 with various types of detectors and various values of other system 
parameters such as dark current and thermal noise in the integrators. 
To apply the Chernov bounds we will need the SIMGF's of the variables 
X at the outputs of the receivers. (See Figs. 3 and 4.) 

An important result needed here is the following: 

Lemma: If C is integer-valued and greater than or equal to zero; and if 
U = ^o Qi where the g { are independent, identically distributed random 
variables (that is, each "count" produced by the C process independently 
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generates g { contributions to the U process through some gain mechanism); 
then the SIMGF of U is 1 

Ms) = fc(*.(s)) (22) 

where \fs (s) is the SIMGF of the identically distributed random variables g, . 

One can model an avalanche photodiode as a photon counter followed 
by a random avalanche multiplier. Each photon counter output "count" 
produces a random number of "counts" at the multiplier output. Thus, 
since we know from Section V the SIMGF of a random multiplier which 
corresponds to \p a (s) in (22), and since we seek the SIMGF of the photo- 
diode output which corresponds to ^pa(s) in (22), it follows that we need 
the SIMGF of the number of counts emitted by a photon counter with 
light incident upon it, which corresponds to ^ c (s) in (22). 

6.4.2 Photon Counter Statistics 

6.4.2.1 Coherent Light. If the light incident upon the photon counter 
is coherent, then it is well known that the SIMGF of the total counts 
emitted in an interval T is given by 6 (see Appendix B) 

tr(s) = [LAM SIC + LAMD] [e* - 1] (23) 

where 

LAMSIG = Total incident light energy i\f Ml 
77 = Detector quantum efficiency 
Ml = Energy of photon at optical frequency used 
LAMD = Mean number of dark current counts before avalanche 
gain in an interval T. 

6.4.2.2 Incoherent Light. If the incident light is incoherent with H 
independent spatial-temporal "degrees of freedom," then the SIMGF of 
the total counts emitted by the photon counter in interval T is (see 
Appendix B) 

\-MMI0 ( ,_ x) ]-»] (24) 



yp c = LAMD[e" - 1] + In 



where LAMSIG is the average total incident energy times y/htt. 
Clearly (24) is the same as (23) as H approaches infinity, which is a 
well known result. 

6.4.3 Final Calculations 

6.4.3.1 Twin-Channel System. Since, for the twin-channel receiver, 
X consists of the difference of two random integrator outputs, we need 
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the following well known result. 4 If X = X x — X 2 , and if Xi and X 2 
are independent, then the SIMGF of X is 

Ms) = iM«) + +xXs). (25) 

Each integrator output contains the sum of the counts emitted by its 
detector and the integral of its thermal noise. The SIMGF of the random 
variable N obtained when Gaussian thermal noise spectral height N Q is 
integrated over an interval T is well known to be 

Ms)=^N T. (26) 

Using (22) through (26) we obtain for the SIMGF of the twin-channel 
receiver output X when the information is in state "one", and the 
optical source is coherent, 



fM ^.V„7' + \^j^ + LAMD 



[M„(s) - 1] 



[yT-Exr rj , , 11; .,;|, 



I I " j^" " + LAMD \[M.(-a) - 1] (27) 



where 



LAMD = mean number of dark current counts before avalanche 
gain in an interval T 
N = spectral height of the thermal noises referred to the 
integrator inputs. 

M„(s) is obtained from (16) through (20) depending upon the particular 
gain mechanism. If the optical source is incoherent, we have 

*x(«) = t'N„T + LAMD[M,(s) + M,(-») - 2] 

We seek the probability that when the information is in state "one", 
X is less than zero, and we therefore decide that the information state 
was "zero." That is, we seek the error probability. One can use (27) 
or (28) and the Chernov bound of (21) to determine the required value 
of LAMSIG = p-T'7i/(hSi) to achieve a desired error probability. 
Since the twin-channel receiver is symmetric, the error probability 
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when the information is in state "zero" is the same as when it is in 
state "one." 

6.4.3.2 Single-Channel System. For the single-channel receiver, we 
need the SIMGF of X under both information states. Call X t the random 
variable X when the information is in the state "one." Call X the 
random variable X when the information is in state zero. Using results 
of Section 6.4, one obtains for coherent light 

4>xM - ^ + \lamd + 2 J^ L ][MM - 1] 



2 



fe.tt - ^ + [lamd + **%**] iMM - i]. 



(29) 



For incoherent light 
2 



<M«) - ^7^ + LAMD[M g (s) - 1] 



+ ln H 1 " P -mi [M ° (S) - 1] T1 



*x.(«) - ^ + LAMD[M g (s) - 1] 



+ ln 



^-^ro^^-iin (30 



One can then use the results of (29) and (30) along with the Chernov 
bounds of (21) to simultaneously find values of LAMSIG = p -T-rj/ihO) 
and the threshold y (see Fig. 4) to ensure some desired error probability 
(which for convenience here will be the same for either information state). 

6.5 Numerical Results 

The Chernov bounds described above were evaluated numerically. 
The results are displayed on the attached figures described below. The 
range of parameter values is realistic and practical, to the best of this 
author's knowledge. The curves presented are those deemed most 
interesting by the author. Other calculations can of course be made. 
Parameters used are defined as follows:* 

LAMSIG = Required mean number of detected photons per 
pulse in the "on" channel of the twin-channel 



* SIG, EXT, G, K, H, and LAMD are input parameters to the program which 
calculates LA MSIG for a desired error rate. 
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receiver or in the "one" state of the single-channel 
receiver. 
LAM SIC EXT = Mean number of detected counts per pulse in the 
"off" channel of the twin-channel receiver or in 
the "zero" state of the single-channel receiver. 
SIG = Normalized thermal noise standard deviation 

= {4/:er/[Re 2 ]}* = {4kdC/e 2 }' 

where e = electron charge, kd = Boltzmann's 
constant • absolute temperature, R = equivalent 
noise resistance at integrator input, T = pulse 
duration, C = T/R = integrator equivalent input 
capacitance. For the results to follow, a reasonable 
value of SIG was chosen to be 6000. 
G = Mean avalanche gain. 
H = Temporal-spatial diversity for incoherent carrier 

case, 
fc = Ratio of ionization probability per unit length of 
weaker and stronger ionizing carriers.* 
LAMD = Dark current counts per interval T before 
avalanche gain. 

Fig. 5 

LAM SIG vs G is plotted for the twin-channel case with k as param- 
eter. SIG was set at 6000, the error rate is 10 _u , LAMD was set at 
5 counts and EXT = 0.01. H = 10,000 which is equivalent to assuming 
a coherent carrier. 

Fig. 6 

The value at optimal gain of LAM SIG vs k is plotted. Points are 
tagged with the optimal G. The receiver is a twin-channel system with 
SIG = 6000, EXT = 0.01, LAMD = 5. H is 10,000 which is equivalent 
to assuming a coherent carrier. The error rate is 10~ B . 

Fig. 7 

LAM SIG vs G is plotted for two values of error rate 10~ 9 and 10" 8 for 



* For these calculations it was assumed that the detector is designed so that the 
stronger ionizing carriers generated optically or associated with dark current enter 
the high-field region from a drift region outside the high-field region. This corresponds 
to initial pairs entering the gain mechanism of x = or x = W as discussed in Section 
III. 
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Fig. 5 — LAMSIG versus gain. 
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Fig. 7 — LAMSIG versus gain. 



a twin-channel system with SIG = 6000, EXT = 0.01, H = 10,000 
LAMD = 5, k = 0.1. 

Fig. 8 

Single- and twin-channel systems are compared. LAMSIG vs G is 
plotted for SIG = 6000, EXT = 0.01, LAMD = 5, H = 10,000, error 
rate = 10~ 9 , fc = 0. Note that from an average power viewpoint the 
single-channel system is 3 dB better than shown if the binary informa- 
tion source is random, since LAMSIG is the energy in "one" state. 

Fig. 9 
Same as Fig. 8 except k = 1. Note the scale change. 

Fig. 10 

LAMSIG vs G for H = 100 and H = 10,000 for twin-channel system. 
SIG = 6000, EXT = 0.01, LAMD = 5, k = 0, error rate = 10" 9 . 

6.6 Further Comments 

When systems were investigated for sensitivity to the choice 
LAMD = 5, EXT = 0.01, it was found that insignificant changes in 
LAMSIG vs G occurred when various combinations of LAMD = 5 
or 50, EXT = 0.01 or 0.001 were tried. 
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Fig. 8 — LAMSIG versus gain. 
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VII. CONCLUSIONS ON APPLICATIONS 



If one assumes that the Chernov bounds are sufficiently tight so that 
actual energy required per bit to achieve specified error rates can be 
compared for various system parameters by comparing the bounds,* 
then one can conclude the following. 



(i) Define k as the ratio of the collision ionization probabilities per 
unit length of the weaker-ionizing to the stronger-ionizing carrier 
(carriers are of course holes and electrons). Assume that the detector 
is designed "well" such that optically and thermally generated carriers 
enter the high-field region from a drift region outside. From the bounds, 
one obtains the result that the required energy per pulse to achieve a 
desired error rate decreases as k decreases for fixed average avalanche 
gain. A value k = is best; but a value k = 0.1 will allow one to operate 
with energy within 10 dB of that required at very high gains with a 
k = device. For each value of k except zero there is an optimal gain 
resulting in minimum required energy per pulse. The optimal gain is 
larger for smaller k. At k = 0.1, the optimal gain is about 100. At 
k = 0, the optimal gain is infinite, but a gain of a few thousand allows 



* For simple cases where both the bounds and actual energy requirements can 
be obtained (for instance for the k = case) the two results differ by a few dB or 
less. 
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close to optimal required energy per pulse. One can conclude that a 
silicon device with k = 0.1 and a gain of about 100 would be a good 
choice for an optical detector. This is true since a detector with k less 
than 0. 1 and yet having gain significantly higher than 100 is not avail- 
able at this time. 

(ii) The required energy per pulse for systems using incoherent 
optical sources differs from that for systems using coherent sources by 
less than a few dB provided the product of the source bandwidth and 
the pulse duration exceeds 100. This is true even if there is no spatial 
incoherence of the light at the detector. 

(Hi) For reasonable parameter values, and assuming a random 
information stream, the single-channel receiver requires about 1.5 dB 
less energy per pulse to achieve a desired error rate than the twin-channel 
receiver. 

(iv) The required energy per pulse is insensitive to reasonable values 
of dark current and extinction ratios. 

(y) For a particular system, a change in the desired error rate from 
10~' J to 10~ 5 results in a change in the required energy per pulse of 
1 to 3 dB, depending upon the avalanche gain. This shows that the 
required energy per pulse is fairly insensitive to the error rate. On the 
other hand, this means that poor error rates will result if insufficient 
loss margin is provided. That is, a small lowering of the received energy 
can greatly increase the error rate. 

APPENDIX A 

In an unpublished work, Mclntyre conjectures (from special case 
calculations) that the probability density of the random gain, defined 
here as p„(n, 0) is given by 

Til n _i_ i L-J/„-*j „-*\"-i 

, m T \r^k + l r (e _ ° 

vM, o) = 



"'^j— g + 2-a 

where k and 8 are the same as in (10) through (13). 

If one makes the assumption that the conjectured p a (n, 0) has sum 
over n normalized to unity for each value of k and for each 5, then one 
obtains the result of (14) by using the definition of the moment-generat- 
ing function and the normalization property. 
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APPENDIX B 

If light of known intensity falls upon a photon counter during an 
interval T, then the probability density of the total number of counts 
emitted is well known 6 to be Poisson distributed as follows 

-(tXMD + A) 

p c (n) = [A + LAMDT - (31) 

Where A* is the total energy incident in the interval T times Ti/hQ,, 
LAMD is the mean number of dark current counts per second times 
the interval T, and r]/hil is the detector quantum efficiency divided by 
the energy in a photon. 

The moment-generating function of the distribution of (31) is given by 

M c (s) = exp [[A + LAMD][e - 1]]. (32) 

If the incident light is a stochastic process, then the moment-generat- 
ing function of the output count distribution is obtained by averaging 
(32) over the probability density of the stochastic total energy incident 
in the interval T 



M 



7 (s) = [ exp ([A + LAMD][e" - l])p(A) dA. (33) 

Jo 



An incoherent light field is normally taken to mean that the complex 
envelope of the classical field is a complex Gaussian random process. 
That is, such a field incident on the photon counter plane can be written 
as 

E( P ,t) = V2re{e(p,t)S ut } 

p t counter plane (34) 

t t (0, T) 

where e(p, t) is a complex Gaussian random process. 

If one expands e(p, t) in its Karhunen-Loeve eigenf unctions, 5 one 
obtains 

«(p, = £ e t t (p, t) 

p e counter plane (35) 

t e (0, T) 



In the text, A is called LAMSIG. 
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where 

f I Mp, m(p,t)d 2 P dt = s*.,- 

''counter plnne •'O 

and the coefficients e k are independent complex Gaussian random vari- 
ables satisfying 

te) = o. 

The energy incident upon the photon counter is 

— A = [ e(p, t)e*(p, I) d 2 p dt= Z \e k \\ (37) 

If one assumes an equal distribution of average energy in roughly H 
"modes," 

7* = 7, 1 ^ k 5g H (3g) 

= 0, k > H 

then it follows that from (33) and the complex Gaussian statistics of 
the e k that 



M c (s) = exp [LAMD(e* - 1)] 



l -v^'-»] H - 



(39) 



Assumption (38) implies that the energy of the incoherent light is 
roughly equally distributed in H degrees of freedom. 
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Previous results on the statistics of avalanche detectors are generalized to 
the case where electrons and holes suffer collision ionizations with unequal 
-probability. It is assumed here that the ratio of collision ionization ■prob- 
abilities per unit length of weaker-to-stronger carrier is a constant k in- 
dependent of position in the high-field region. The moment-generating 
function of the random avalanche gain G is obtained as a function of k 
and the average gain G, and is used to obtain Chernov bounds on error rates 
of digital optical receivers employing avalanche detectors. It is shown the 
required energy per pulse to achieve a given error rate decreases as k de- 
creases for fixed G. For each k > 0, there is an optimal mean gain G op , 
resulting in minimum required energy per pulse. At k = 0.1, G opl tt 100 
and the required energy is within 10 dB of that required with very high 
gains (a few thousand) at k = 0. 

I. INTRODUCTION 

In a previous paper 1 results on the statistics of two particular ava- 
lanche detectors with applications to optical communication were 

3075 
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presented. It was required either that only one carrier suffer collision 
ionizations in the high-field region (unilateral gain) or that both carriers 
suffer collision ionizations with equal probability per unit length in the 
high-field region. The present work allows for more general unequal 
ionization probabilities per unit length with the requirement only that 
the ratio of the two quantities be constant throughout the high-field 
region. The moment-generating function of the random gain is obtained 
as a function of this ratio and the average gain. The results are con- 
sistent with unpublished conjectures of R. J. Mclntyre. 2 The moment- 
generating functions are used to obtain Chernov bounds on the error 
rates of digital optical receivers employing avalanche detectors and 
using either coherent or incoherent light. Results on avalanche statistics 
are summarized in Section V. Numerical results on the Chernov bounds 
are given in Section VI (6.5) and Section VII. 

II. MODEL OF THE AVALANCHE DETECTOR 

The avalanche detector is a device in which thermally or optically 
generated hole-electron pairs generate additional hole-electron pairs 
through collision ionizations. Within the device there is a "high-field 
region" where holes have probability j3(x) per unit length (which depends 
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Fig. 1 — Avalanche detector. 
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upon the position x) of suffering a collision ionization as they travel to 
the left under the influence of the electron field (see Fig. 1). Electrons 
traveling to the right have a probability a(x) per unit length of collision 
ionization. Carriers can be created within the high-field region due to 
thermal effects or due to the presence of incident light. Carriers can 
also drift into the region if they are generated outside of the region. It is 
assumed that all collision ionizations are independent. This requires 
that the mean distance between ionizing collisions be large compared 
to the distance over which a carrier can randomize its momentum after 
a collision. 3 Hole-electron pairs created through collision ionization can 
in turn generate additional pairs by the same mechanism. This is the 
avalanche process. 

III. THE STATISTICS 

We seek the statistics of the random total number of hole-electron 
pairs which result ultimately through collision ionizations when an 
initial hole-electron pair is injected into the high field at some position x. 
Define p a (n, x) as the probability that n pairs ultimately result in- 
cluding the initially injected pair. The moment-generating function of 
the number of pairs M„(s) is therefore 4 

M„(s, x) = £ VM, x)e n . (1) 

i 

We shall derive M g (s, x). Before proceeding we must review some well 
known results which will be needed in that derivation. 

If { Xi \ are random variables which are independent, then the moment- 
generating function of the sum of the {a;,} is the product of the individual 
moment-generating functions. 4 The semi-invariant moment-generating 
function SIMGF of a random variable X having probability density 
p x (x) is the natural logarithm of the moment-generating function of X 



Ms) - In [M x (s)] = In [ f) p*(*)e M ] • 



(2) 



The SIMGF of a sum of n independent random variables is the sum of 
the individual SIMGF's. 

We can now proceed to derive M„(s, x). Divide the high-field region 
into K intervals of width dX = W/K. See Fig. 2. Label these intervals 
1, 2, 3, • • • j, • • • K. If a hole-electron pair is injected into interval j, 
define, as above, the probability density of the total number of pairs 
ultimately resulting in the avalanche process (including the initial pair) 
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as p (n, x) where x is taken as the center of interval j. The hole of the 
initial pair moves to the left and the electron to the right toward 
x = or x = W respectively. As they pass through their respective 
intervals, new pairs may be created in each through collision ionization. 
We shall assume that the interval width dX is sufficiently narrow so 
that the initial pair carriers create either one or no new pairs in each 
interval. If the initial hole or electron generates a new pair in some 
interval fc, then that new pair will ultimately generate N k pairs including 
itself through the avalanche process. Thus with each interval we can 
associate a number of pairs N k . This number equals zero if the appro- 
priate initial pair carrier suffers no collision ionizations in interval k. 
This number equals one or more if the appropriate initial pair carrier 
suffers a collision ionization in interval k. The total number of pairs 
ultimately generated through the avalanche process including the 
initial pair is one plus the sum of the {N k \. Since collision ionizations 
are all independent, all the N k are independent. Thus we have the 
SIMGF of the total number of pairs given by 



^(s, x) = s + 2 *N t (s) 



(3) 



where s is the SIMGF of the deterministic initial pair and ^ , A r l (s) is 
the SIMGF of N k . 

The SIMGF of N k is obtained as follows. The probability that 
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N k = is 1 — 7 dX where 7 = a(x) if interval k is to the right of interval 
;' where the initial pair enters, y = /3(x) if k is to the left of ;. (x is eval- 
uated at the center of interval k.) The probability that N k = Z > is 
the probability of a collision ionization in interval k times the probability 
that a new pair created at interval k ultimately results in Z pairs 
including itself. That is 



Pr (N k = Z) = 7 dX p g (Z, x) for Z > 
with x evaluated at the center of interval k. Thus 



(4) 



l_z=o 



lM«) = ln LPr (Nk = Z)e 



= In 



(i-7 dxy° + 2 y dx Pb (z, x y 



(5) 



where 



= In [1 - 7 dX + 7 dX M„(s, x)] 
= ln[l - t^ + 7^^"' 1 '] 

7 = <x(z) if k > j 

7 = /3(z) if k < j 

x = center of interval k. 

Using (3) and (5) and taking the limit as dX gets infinitely small, one 
obtains 



h(s,x) =s+ / /3(x')[c*' (, - z,) - 1]^' 

Jo 

+ J a(x')(e*'"""' - 1) dx'. (6) 

Equation (6) is the critical equation for determining ^„(s, x), and thus 
il/ s (s, x) = exp [^„(s, a;)]. Using Leibnitz's rule for differentiation of 
integrals one obtains 



dx 



*.(«, X) = \fi(x) - a(xW> 



~ 1). 



The solution of (7) is 
f Q (s, x) = In 



1 



1 - C exp I / Q5(x') - a(x')] dx') 



(7) 



(8) 
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where C = (e*' (,,0) - l)/e*» <,,0) which can be checked by substitution. 
Substituting (8) into (6) one obtains the particular result 






f «(*') 

•'o 



- 1 



1 - C exp / \0{x") - a(x")] dx" 



dx'. (9) 



If one makes the assumption that at each point in the high-field 
region 

j8(x) = k-a(x) (10) 

where k is a constant, one can solve (9) to obtain 

*.(«, 0) = s - 8 + — - In [- 



* - 1 Ulf„(«, 0) - e (t - 1)5 [M e (s, 0) - 1]J 

= 8 + j4-fc ln f^-(«» °) - e ( *" ,)S [^(s, 0) - 1]] (11) 



where 

M, 



(s, 0) = e*' l '- 0) and S = f a(x) dx. 

Jo 



One can write (11) in another way by making a substitution.* Define 
8'(s) implicitly by 

e 5 ' (8) = M„(s,0)e s - y (12) 

that is, 

e #.c..o) = Mb(s>0 ) = e- ,+l ' ( '». (13) 

Using (12) in (11) one obtains the implicit equation 



-mm.) _ e -vM = e » [g -*« _ e -ij (14) 



which determines 8'(s) and thus M B (s, 0) through (13). 

Equation (14) is still not explicit. A numerical technique for solution 
is discussed in the next section. One can use (11) [or (14)] and (8) to 
obtain $ g (s, x) or M g (s, x) for any x. Recall that x is the point of entry 
of the initial pair. In the applications we shall be concerned with x = 
or x = W. That is, pairs are generated in a drift region outside the high- 
field region with carriers drifting into the high-field region. 



* Equation (14) will follow from (11) and (12) by tedious algebra. Further, (14) 
will not be used in the following results except to compare with Mclntyre's work 
in Appendix A. 
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The equation (14) is consistent with some unpublished conjectures of 
Mclntyre given in Appendix A. 

IV. NUMERICAL SOLUTION'S 

Equations (11) through (14) can be solved numerically. One technique 
is to differentiate (11) to obtain the result 

f a M,(s, 0) = M B (s, 0)(^-=^)[l - \ (M e (s, 0)<fVr '] (15) 

where M o (0, 0) = 1. 

Equation (15) can be integrated with a computer to obtain M„(s, 0) 
explicitly. 

V. SUMMARY OF ANALYTIC RESULTS ON AVALANCHE STATISTICS 

From Sections I through IV and a previous paper by this author, 1 
we obtained the following: 

Assumptions: 

Holes travel toward x = 0, electrons toward x = W. 
Hole ionization probability per unit length = fi(x). 
Electron ionization probability per unit length = <x(x). 
0(x) = k-a(x), k a constant for all x. 
High-field region width = W. 

Definitions: 

f W 

8 = I a(x) dx, 

Jo 

p„(n, x) = probability that if an initial pair enters the high-field region 
at point x, n pairs will ultimately result through the avalanche process 
including the initial pair, 

00 

G(x) = mean avalanche gain = 23 n Vo( n > ■ r )» 

i 

M g (s, x) = moment-generating function of p„(n, x) 

n-l 

Results: 
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L M - (s - * - 1 1 - (M £(!lo7 1} exp [ (A - :) f aM "IT 

for all fc. (16) 
For k = (Unilateral Gain) 

2a. Af.(«, 0) = [1 - e 6 [l - e - ']] -1 where (7(0) - e\ (17) 

2b. p,(n, 0) = | (^ 1 )" where G = (7(0). (18) 

For k ^ 0, k ^ 1 

3. ^Af,(s, 0) = M,(s, 0)(^ I )[l - | (ilf,(«, O^-V)*" 1 ] ', 

where 

M,(0, 0) = 1 (19) 

where 

^>-(^)[ 1 -r" , -"j"' 

For A; = 1 (Equal Ionization) 

4. ^Mi( S , 0) = M.(s, 0)[l - (^=^)ilf,(s, 0)] 

M B (0, 0) = 1 
where 

G = (7(0). (20) 

VI. APPLICATIONS TO RECEIVERS USING AVALANCHE DETECTORS 

6.1 General Comments 

We shall next apply the results of Section V to obtain bounds on the 
error rates of digital receivers. The receivers to be discussed here are 
the single- and twin-channel systems described below. We shall upper- 
bound the power required at the receiver to obtain a desired error rate 
using the Chernov bounds. 

6.2 The Receivers 

The twin-channel receiver is shown in Fig. 3. Depending on the state 
of a binary information source, one of two channels has optical output 
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Fig. 3 — Twin-channel system. 



power p for duration T while the other has output power p-EXT for 
duration T. EXT signifies extinction. EXT would ideally be zero, but 
is left finite and less than unity for practical reasons. The optical power 
falling on an avalanche photo diode causes the emission of photo- 
electrons which are multiplied along with the detector "dark current" 
through the avalanche gain mechanism. The detector outputs (multi- 
plied counts) are integrated in devices having thermal noises referred to 
their respective inputs. The integrator outputs are subtracted and the 
difference X is compared to a threshold of zero to decide what the 
information state was. 

The single-channel system is essentially half the twin-channel system 
as shown in Fig. 4. The single integrator output X is compared to a 
threshold y to decide upon the information state. 

6.3 The Chernov Bounds 

The Chernov bound is a useful tool for bounding the probability that 
a random variable X will lie above or below a given threshold 7. It is 
given as follows 5 

Pr(X> 7 ) £«"<•»-♦*<•> | w ; (1) 

provided s > 0, 
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Pr (X <y)^ e *'M-'**M | t-#;w 

provided s < 0, (21) 



where 



and 



Ms) = SIMGF of X m In \f p x (x)e" dx\ 



rx{s) = y 8 m*)- 



6.4 Chernov Bounds for the Two Systems 

6.4.1 Preliminaries 

We wish to determine the required power p to achieve a desired 
decision error probability for each of the receivers discussed in Section 
6.2 with various types of detectors and various values of other system 
parameters such as dark current and thermal noise in the integrators. 
To apply the Chernov bounds we will need the SIMGF's of the variables 
X at the outputs of the receivers. (See Figs. 3 and 4.) 

An important result needed here is the following: 

Lemma: If C is integer-valued and greater than or equal to zero; and if 
U = 2^o 9< where the g { are independent, identically distributed random 
variables {that is, each "count" produced by the C process independently 
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generates g t contributions to the U process through some gain mechanism); 
then the SIMGF of U is 1 

Ms) = iMlMs)) (22) 

where ^„(s) is the SIMGF of the identically distributed random variables g, . 

One can model an avalanche photodiode as a photon counter followed 
by a random avalanche multiplier. Each photon counter output "count" 
produces a random number of "counts" at the multiplier output. Thus, 
since we know from Section V the SIMGF of a random multiplier which 
corresponds to i^ 3 (s) in (22), and since we seek the SIMGF of the photo- 
diode output which corresponds to \pu (s) in (22), it follows that we need 
the SIMGF of the number of counts emitted by a photon counter with 
light incident upon it, which corresponds to \f/ c (s) in (22). 

6.4.2 Photon Counter Statistics 

6.4.2.1 Coherent Light. If the light incident upon the photon counter 
is coherent, then it is well known that the SIMGF of the total counts 
emitted in an interval T is given by 6 (see Appendix B) 

Ms) = [LAM SIC + LAMD] [e' - 1] (23) 

where 

LAMSIG = Total incident light energy ij/Mt 
7} = Detector quantum efficiency 
Ml = Energy of photon at optical frequency used 
LAMD = Mean number of dark current counts before avalanche 
gain in an interval T. 

6.4.2.2 Incoherent Light. If the incident light is incoherent with H 
independent spatial-temporal "degrees of freedom," then the SIMGF of 
the total counts emitted by the photon counter in interval T is (see 
Appendix B) 



tc = LAMD[c" - 1] + In 



, LAMSIG . . 1N 

1 - jj (e - 1) 



(24) 



where LAMSIG is the average total incident energy times r?/M2. 
Clearly (24) is the same a.s (23) as H approaches infinity, which is a 
well known result. 

6.4.3 Final Calculations 

6.4.3.1 Twin-Channel System. Since, for the twin-channel receiver, 
X consists of the difference of two random integrator outputs, we need 
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the following well known result. 4 If X = X, — X 2 , and if X x and X 2 
are independent, then the SIMGF of X is 

Ms) = *x,(«) + *r.(-«). (25) 

Each integrator output contains the sum of the counts emitted by its 
detector and the integral of its thermal noise. The SIMGF of the random 
variable N obtained when Gaussian thermal noise spectral height iV is 
integrated over an interval T is well known to be 

+ N (s)= S ^N T. (26) 

Using (22) through (26) we obtain for the SIMGF of the twin-channel 
receiver output X when the information is in state "one", and the 
optical source is coherent, 

Ms) = s 2 N T + fc|^ + LAMDpfM - 1] 



H 



I'' *'%* Tr, + LAMD 



[M,(-s) - 1] (27) 



where 



LAMD = mean number of dark current counts before avalanche 
gain in an interval T 
N = spectral height of the thermal noises referred to the 
integrator inputs. 

M B (s) is obtained from (16) through (20) depending upon the particular 
gain mechanism. If the optical source is incoherent, we have 

Ms) = s 2 N T + LAMD[M B (s) + M,(-s) - 2] 

+ in ib - w ^> - «n 

+ i 4t i - p - z li^^ ( - s) - i] ]" / 1- (28) 

We seek the probability that when the information is in state "one", 
X is less than zero, and we therefore decide that the information state 
was "zero." That is, we seek the error probability. One can use (27) 
or (28) and the Chernov bound of (21) to determine the required value 
of LAMSIG = ■p-T--r)f(hQ) to achieve a desired error probability. 
Since the twin-channel receiver is symmetric, the error probability 
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when the information is in state "zero" is the same as when it is in 
state "one." 

6.4.3.2 Single-Channel System. For the single-channel receiver, we 
need the SIMGF of X under both information states. Call X x the random 
variable X when the information is in the state "one." Call X the 
random variable X when the information is in state zero. Using results 
of Section 6.4, one obtains for coherent light 

**.(*) = ^f^ + \_LAMD + 2 ^ 3 ][M g (8) - 1] 

feW = ^t + [LAMD + p - T '*£ T '* \lM.(8) - 1]. (29) 
For incoherent light 
*x.« = ^f~ + LAMD[M B (s) - 1] 



+in [[ i -w ?[M ' (s) - i] r] 



*x.(«) = ^^ + LAMD[M.(*) ~ 1] 



^[fr-*^™ -«]-]■ 



(30) 



One can then use the results of (29) and (30) along with the Chernov 
bounds of (21) to simultaneously find values of LAMSIG = p-T- ??/(M2) 
and the threshold y (see Fig. 4) to ensure some desired error probability 
(which for convenience here will be the same for either information state). 

6.5 Numerical Results 

The Chernov bounds described above were evaluated numerically. 
The results are displayed on the attached figures described below. The 
range of parameter values is realistic and practical, to the best of this 
author's knowledge. The curves presented are those deemed most 
interesting by the author. Other calculations can of course be made. 
Parameters used are defined as follows:* 

LAMSIG = Required mean number of detected photons per 
pulse in the "on" channel of the twin-channel 



* SIG, EXT, G, K, H, and LAMD are input parameters to the program which 
calculates LA MSIG for a desired error rate. 
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receiver or in the "one" state of the single-channel 
receiver. 
LAM SIG • EXT = Mean number of detected counts per pulse in the 
"off" channel of the twin-channel receiver or in 
the "zero" state of the single-channel receiver. 
SIG = Normalized thermal noise standard deviation 

= {4/:07V[Re 2 ])* = {IkdC/e 2 }* 

where e = electron charge, fid = Boltzmann's 
constant • absolute temperature, R = equivalent 
noise resistance at integrator input, T = pulse 
duration, C = T/R = integrator equivalent input 
capacitance. For the results to follow, a reasonable 
value of SIG was chosen to be 6000. 
G = Mean avalanche gain. 
H = Temporal-spatial diversity for incoherent carrier 

case. 
k = Ratio of ionization probability per unit length of 
weaker and stronger ionizing carriers.* 
LAMD = Dark current counts per interval T before 
avalanche gain. 

Fig. 5 

LAM SIG vs G is plotted for the twin-channel case with k as param- 
eter. SIG was set at 6000, the error rate is 10~ n , LAMD was set at 
5 counts and EXT = 0.01. H = 10,000 which is equivalent to assuming 
a coherent carrier. 

Fig. 6 

The value at optimal gain of LAM SIG vs fc is plotted. Points are 
tagged with the optimal G. The receiver is a twin-channel system with 
SIG = 6000, EXT = 0.01, LAMD = 5. H is 10,000 which is equivalent 
to assuming a coherent carrier. The error rate is 10~ 9 . 

Fig. 7 

LAM SIG vs G is plotted for two values of error rate 10 -9 and 10 -5 for 



* For these calculations it was assumed that the detector is designed so that the 
stronger ionizing carriers generated optically or associated with dark current enter 
the high-field region from a drift region outside the high-field region. This corresponds 
to initial pairs entering the gain mechanism of x = or x = W as discussed in Section 
III. 
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Fig. 5 — LAMSIG versus gain. 
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Fig. 6 — LAMSIG at optimal gain versus k. 
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Fig. 7 — LAMSIG versus gain. 



a twin-channel system with SIG = 6000, EXT = 0.01, H = 10,000 
LAMD m 5, k m 0.1. 

Fig. 8 

Single- and twin-channel systems are compared. LAMSIG vs G is 
plotted for SIG = 6000, EXT = 0.01, LAMD = 5, H = 10,000, error 
rate = 10~ 9 , k = 0. Note that from an average power viewpoint the 
single-channel system is 3 dB better than shown if the binary informa- 
tion source is random, since LAMSIG is the energy in "one" state. 

Fig. 9 
Same as Fig. 8 except k = 1. Note the scale change. 

Fig. 10 

LAMSIG vs G for H = 100 and H = 10,000 for twin-channel system. 
SIG = 6000, EXT = 0.01, LAMD = 5, k = 0, error rate = 10" 9 . 

6.6 Further Comments 

When systems were investigated for sensitivity to the choice 
LAMD = 5, EXT = 0.01, it was found that insignificant changes in 
LAMSIG vs G occurred when various combinations of LAMD = 5 
or 50, EXT = 0.01 or 0.001 were tried. 



STATISTICS OF AVALANCHE DETECTORS 



3091 



3000 r 


1 








2500 








SIG=6000 
LAMD=5 


2000 








EXT = 0.0l 

H = 10 4 (COHERENT) 

k=o 

ERROR RATE=10" 9 


O 










2 1500 
< 

_i 


1 SINGLE 


CHANNEL 




1000 


" 1 ^,,-TWIN 


CHANNEL 




500 


1 











1 1 


1 1 1 -I 



250 500 750 1000 1250 1500 1750 2000 

G 



Fig. 8— LAMSIG versus gain. 
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Fig. 9 — LAMSIG versus gain. 
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Fig. 10 — LAMSIG versus gain. 



VII. CONCLUSIONS ON APPLICATIONS 



If one assumes that the Chernov bounds are sufficiently tight so that 
actual energy required per bit to achieve specified error rates can be 
compared for various system parameters by comparing the bounds,* 
then one can conclude the following. 



(i) Define k as the ratio of the collision ionization probabilities per 
unit length of the weaker-ionizing to the stronger-ionizing carrier 
(carriers are of course holes and electrons). Assume that the detector 
is designed "well" such that optically and thermally generated carriers 
enter the high-field region from a drift region outside. From the bounds, 
one obtains the result that the required energy per pulse to achieve a 
desired error rate decreases as k decreases for fixed average avalanche 
gain. A value k = is best; but a value k = 0.1 will allow one to operate 
with energy within 10 dB of that required at very high gains with a 
k = device. For each value of k except zero there is an optimal gain 
resulting in minimum required energy per pulse. The optimal gain is 
larger for smaller k. At k = 0.1, the optimal gain is about 100. At 
k = 0, the optimal gain is infinite, but a gain of a few thousand allows 



* For simple cases where both the bounds and actual energy requirements can 
be obtained (for instance for the k = case) the two results differ by a few dB or 
less. 
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close to optimal required energy per pulse. One can conclude that a 
silicon device with k = 0.1 and a gain of about 100 would be a good 
choice for an optical detector. This is true since a detector with k less 
than 0. 1 and yet having gain significantly higher than 100 is not avail- 
able at this time. 

(it) The required energy per pulse for systems using incoherent 
optical sources differs from that for systems using coherent sources by 
less than a few dB provided the product of the source bandwidth and 
the pulse duration exceeds 100. This is true even if there is no spatial 
incoherence of the light at the detector. 

(Hi) For reasonable parameter values, and assuming a random 
information stream, the single-channel receiver requires about 1.5 dB 
less energy per pulse to achieve a desired error rate than the twin-channel 
receiver. 

(w) The required energy per pulse is insensitive to reasonable values 
of dark current and extinction ratios. 

(v) For a particular system, a change in the desired error rate from 
10~ 9 to 10~ 5 results in a change in the required energy per pulse of 
1 to 3 dB, depending upon the avalanche gain. This shows that the 
required energy per pulse is fairly insensitive to the error rate. On the 
other hand, this means that poor error rates will result if insufficient 
loss margin is provided. That is, a small lowering of the received energy 
can greatly increase the error rate. 

APPENDIX A 

In an unpublished work, Mclntyre conjectures (from special case 
calculations) that the probability density of the random gain, defined 
here as p a (n, 0) is given by 

r fe + !>-(.- - e -V- 

p B {n, 0) = -. 

where k and 5 are the same as in (10) through (13). 

If one makes the assumption that the conjectured j> g {n, 0) has sum 
over n normalized to unity for each value of k and for each 5, then one 
obtains the result of (14) by using the definition of the moment-generat- 
ing function and the normalization property. 
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APPENDIX B 

If light of known intensity falls upon a photon counter during an 
interval T, then the probability density of the total number of counts 
emitted is well known 6 to be Poisson distributed as follows 

-ILAMD + S.) 

p c (n) = [A + LAMD] n - - (31) 

Where A* is the total energy incident in the interval T times rj/hto, 
LAMD is the mean number of dark current counts per second times 
the interval T, and rj/hQ, is the detector quantum efficiency divided by 
the energy in a photon. 

The moment-generating function of the distribution of (31) is given by 

M c (s) = exp [[A + LAMD][e' - 1]]. (32) 

If the incident light is a stochastic process, then the moment-generat- 
ing function of the output count distribution is obtained by averaging 
(32) over the probability density of the stochastic total energy incident 
in the interval T 



M, 



:(s) = f exp ([A -f- LAMD][e a - l])p(A) dA. (33) 

Jo 



An incoherent light field is normally taken to mean that the complex 
envelope of the classical field is a complex Gaussian random process. 
That is, such a field incident on the photon counter plane can be written 
as 

E( P ,t) = V2rele( P ,t)e iU \ 

p c counter plane (34) 

t e (0, T) 

where e(p, t) is a complex Gaussian random process. 

If one expands e(p, in its Karhunen-Loeve eigenfunctions, 5 one 
obtains 

«(p, = £ e k (p k (p, t) 

p t counter plane (35) 

t t (0, T) 



* In the text, A is called LAMSIG. 
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where 

[ [ <t> k (p,t)< P * j {p,t)d 2 pdt = h.i 

•'counter plane J 

and the coefficients e k are independent complex Gaussian random vari- 
ables satisfying 

(e k e*l) = y k 8 k ,i ,^Q) 

(e k e,) = 0. 
The energy incident upon the photon counter is 

^ A = [ e(p, t)e*(p, t) d 2 p dt= £ \e k \ 2 . (37) 

If one assumes an equal distribution of average energy in roughly H 
"modes," 

7, = 7, 1 ^ k ^ H (3g) 

= 0, k > H 

then it follows that from (33) and the complex Gaussian statistics of 
the e k that 



M c (s) = exp [LAMD{e - l)]-[l - ^7^ - 1)J • 



(39) 



Assumption (38) implies that the energy of the incoherent light is 
roughly equally distributed in H degrees of freedom. 
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